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Holes doped within the square CuC>2 network specific to the cuprate superconducting materials 
have oxygen 2p character. We investigate the basic properties of such oxygen holes by wavefunction- 
based quantum chemical calculations on large embedded clusters. We find that a 2p hole induces 
ferromagnetic correlations among the nearest-neighbor Cu 3d spins. When moving through the anti- 
ferromagnetic background the hole must bring along this spin polarization cloud at nearby Cu sites, 
which gives rise to a substantial reduction of the effective hopping parameters. Such interactions 
can explain the relatively low values inferred for the effective hoppings by fitting the angle-resolved 
photoemission data. The effect of the background antiferromagnetic couplings of renormalizing the 
effective nearest-neighbor hopping is also confirmed by density-matrix renormalization-group model 
Hamiltonian calculations for chains and ladders of CuCU plaquettes. 



I. INTRODUCTION 

Most theoretical models for high-temperature super- 
conductivity are based on Hubbard or t — J Hamilto- 
nians with one orbital per CUO4 plaquette. For the 
square lattice, the tight-binding part of such models reads 
in the k representation e(k) = — 2t(cosk x + cosk y ) + 
it' cos k x cos k y — 2t" (cos2k x + cos2k y ) + ... , where t is 
the hopping integral between nearest neighbors and t' 
and t" are hoppings between second and third order 
neighbors. The values of these parameters were esti- 
mated by fitting the angle-resolved photoemission spec- 
troscopy (ARPES) experimental data, see for example 
Refs. [l|, 0, Q , and also by first principles investigations, 
both periodic density-functional calculations within the 
local density approximation (LDA) 0, H, @ and em- 
bedded cluster, wavefunction-based calculations 0, H, Q . 
The analysis of the LDA conduction bands in several cop- 
per oxide compounds yields a value of 0.4-0.5 eV for t and 
a ratio between the nearest-neighbor and next-nearest- 
neighbor hopping integrals t' /t~0.\5 for La2Cu04 0, @| 
and 0.33 for Tl 2 Ba 2 Cu06 6]. Wavefunction-based, quan- 
tum chemical calculations on finite clusters predict very 
similar nearest-neighbor hoppings, with t varying be- 
tween 0.5 and 0.6 eV for different cuprate superconduc- 
tors @, H, @ , with the observation that t correponds in 
this case to the propagation of an oxygen 2p hole and 
not to d-like conduction-band states. However, these es- 
timates are quite far from the numbers obtained by fitting 
the ARPES data. Such procedures give usually nearest- 
neighbor hopping integrals in the range of 0.10-0.30 eV 
for hole doping and 0.10-0.20eV for the electron doped 
materials [2,0>O|. 

This paper presents a careful analysis of the character- 
istics of oxygen 2p holes doped into the spin- 1/2 antifer- 
romagnetic Cu02 layer [1(1. W e perform ab initio multi- 
configuration calculations [ljj on clusters that are large 



enough to account for charge and spin relaxation and po- 
larization effects in the near surrounding. We show that 
the hopping of the doped particle is strongly renormal- 
ized by spin interactions with and among the 5 = 1/2 Cu 
of neighbors and that by properly taking into account 
these effects, good agreement is found with the estimates 
extracted from the photoemission data for the hopping 
integrals. Our study is able to explain thus the difference 
between these "fitted" tight-binding parameters and the 
results of previous embedded cluster quantum chemical 
investigations or periodic density- functional calculations. 
The strong renormalization of the nearest-neighbor effec- 
tive hopping found by ab initio wavefunction-based quan- 
tum chemical methods is confirme d by density-matrix 
renormalization-group (DMRG) EH [Of model Hamilto- 
nian calculations on chains and ladders of CUO4 plaque- 
ttes. 



II. CHARGE AND SPIN CONFIGURATION 
AROUND A DOPED HOLE 

To illustrate the nature of an oxygen hole in the CuC>2 
plane, we first discuss results of multiconfiguration self 
consistent field (MCSCF) calculations on a 9-plaquette 
square cluster, see Fig.l. This cluster includes nine Cu 
ions and twenty four in-plane oxygens but no apex lig- 
ands. The geometrical structure corresponds to that of 
the hole doped Lai.gsSri.gsCuC^ compound [l4j]. We 
use effective core potentials (ECPs) for the Is, ...3s inner 
shells of each transition metal center and the Is core of 
the oxygens at the boundaries of the cluster, as developed 
by Seijo et al. [H| and Bergner et al. [HI, respectively. 
For the outer shells of these ions we apply the following 
gaussian-type basis sets: Cu (9s6p6d) /[3s3p3d] [15| and 
O (4s5p)/[2s3p] [Hj]. Still, the rest of the ligands, i.e. 
those bridging the cluster Cu ions (see Fig. 1(a)), are rep- 
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FIG. 1: (a) The 9-plaquette [CU9O24] cluster employed for 
studying the distribution of an 2p hole. O ions that are 
modeled with all-electron basis sets are shown in grey, (b)-(c) 
B and AB pd orbitals defining the ZR-like state on the central 
Cu c 04 plaquette. Their occupation numbers are 1.80 and 
0.20. Positive and negative lobes of the p and d functions are 
shown in white and gray, respectively. For the B combination 
the p and d AOs are strongly overlapping. 



resented with all-electron basis sets with a (14s9p)/[4s3p] 
contraction scheme [l7| . Our clusters are always embed- 
ded in large arrays of point charges at the experimen- 
tal lattice positions that reproduce the crystal Madclung 
field. The Cu 2+ and La 3+ nearby neighbors are repre- 
sented by total ion potentials • The calculations were 
performed with the MOLCAS 6 program package [l9[ . 

Multiconfiguration, complete active space (CAS) [ll| 
calculations were carried out for a single hole "doped" 
into the 9-plaquette cluster. We employed a minimal 
CAS for constructing the multiconfiguration wavefunc- 
tion, with ten electrons distributed in all possible ways 
over ten active orbitals. Those are the nine Cu 3d x 2_ y 2 
orbitals plus one O 2p orbital. The rest of the Cu 3d and 
O 2p orbitals and the other lower-energy orbitals are dou- 
bly occupied in all configurations and do not participate 
to the construction of the MC expansion. The ground- 
state has overall singlet spin multiplicity and corresponds 
to the formation of a pd Zhang- Rice (ZR) 20] -like config- 
uration on the central Cu c 04 plaquette. Mulliken charge 
and spin populations, MCPs and MSPs, respectively, are 
collected for this state in Table I. The bonding (B) and 
antibonding (AB) pd ZR orbitals are plotted in Fig.l(b- 
c). Their occupation numbers are 1.80 and 0.20, respec- 
tively. The Mulliken charges of the oxygen 2p x and 2p y 
atomic orbitals (AOs) forming the pd B and AB com- 
binations on the ZR plaquette are similar to the values 
reported in [2l], [22J. The Mulliken charges associated 
with the Cu 3d x 2_ y 2 functions are slightly larger due to 
the different basis sets employed here. 

The formation of CUO4 pd singlet states as predicted 
by Zhang and Rice [20] ] was previously confirmed by the 
ab initio study of Calzado and Malrieu 0, Q ■ Two of us 



investigated in Refs. [2l|, |22| electron-lattice interactions 
associated with such configurations. It was found that 
the ZR state is stabilized by lattice relaxation effects in- 
volving a shortening by few percents of the Cu-0 bonds. 
The magnitude of the stabilization energy indicates that 
the doped holes form small polarons in cuprates, at least 
at low concentrations. In this paper, however, we are 
mainly concerned with pure electronic effects, that deter- 
mine the characteristics of the ZR-like quasiparticle. All 
calculations are therefore performed on high-symmetry, 
undistorted clusters using the average Cu-0 distances 
reported in Ref. [lij ]. For the [Cug024] square cluster de- 
picted in Fig.l, both ZR-like states and broken-symmetry 
solutions with the 2p hole having the largest weight at a 
single oxygen site were obtained. Such O 2p 5 -like states 
were discussed in Ref. [21]. The energy difference be- 
tween the ZR Cu c 04 state and a broken-symmetry O c 
2p 5 -like state is 5 meV in the [CU9O24] cluster, with the 
ZR configuration lower in energy. The broken-symmetry 
configurations may be favored, however, by inter-carrier 
Coulomb and spin interactions. Effects of such inter- 
actions between two or more holes were investigated in 
Ref. [23| by calculations on 7-plaquette linear clusters. 
For an undistorted 7-plaquette cluster and a single doped 
hole, the energy difference between the ZR state and a 
broken-symmetry 2p 5 -like state, both involving oxygens 
ions on the same plaquette, is again very small, about 5 
meV, but the order of the states is reversed. This differ- 
ent ordering of the states is due to finite-size effects in 
the linear cluster [24| • To conclude this part of our dis- 
cussion, the calculations on the large 9-plaquette cluster 
show that an isolated O 2p hole, i.e. an oxygen hole at 
sufficiently low doping, should be viewed as part of a pd 
ZR-like state even if the local lattice relaxation effects 
are neglected. Our estimation of the nearest-neighbor 
and next-nearest-neighbor hopping matrix elements will 
be based on such a ZR-like quasiparticle picture. 

A distinctive feature of the results listed in Table I is 
the nature of some of the spin couplings. The spin pop- 
ulations given in the last column are obtained from the 
difference between the "up" and "down" spin densities. 
The very low values associated with the Cu d x 2_ y 2 and 
O p x and p y orbitals forming the linear combinations de- 



TABLE I: Mulliken population analysis illustrating the dis- 
tribution of an O 2p hole and the nature of the Cu-Cu spin 
couplings, see text. CASSCF results for a 9-plaquette cluster. 
Notations as in Fig.l are used. For each of the O c , Cu nn , and 
Outran ions, there are another three equivalent sites. 



Relevant AOs 


MCPs 


MSPs 


3^3-2 _y2 


1.17 


0.06 




1.62" 


-0.01 


CU n n Sd x 2_y2 


1.27 


0.31 


CU^J^ Sdj,2_y2 


1.28 


-0.32 



"The Mulliken charges of the other in-plane oxygen atomic or- 
bitals arc not smaller than 1.8. 
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TABLE II: Nearest-neighbor hopping integrals as obtained by 
CASSCF/CASSI calculations on clusters of different sizes, see 
text. Mulliken charge populations for the O 2p AOs forming 
the a p x ,y-d x 2_ y 2 bonds on a ZR plaquette are also listed. 



Cluster 


t(eV) 


O-hole AOs, MCPs: 
Pv, Py, P* 


[CueOis]" 


0.310 


1.61, 1.61, 1.60 


[CueOir] 6 


0.165 


1.61, 1.63, 1.63 




0.160 


1.60, 1.64, 1.63 


[Cuio0 2 9] c 


0.135 


1.65, 1.61, 1.63 



"Linear cluster, all-electron (AE) basis sets (BSs) 28] . 
b 3 by 2 plaquettes, see text. Results with AE BSs, first line, and 
with ECPs for the ions outside the active region, next line. 
c ECPs were used for part of the Cu and O ions, see text. 



picted in Fig.l(b-c) indicate indeed that the dominant 
contribution to the many-electron wavefunction comes 
from configurations where the d and p holes on the Cu c 04 
plaquette are coupled to a singlet. The fact that these 
numbers are not exactly zero shows that configurations 
related to different and more complicated spin coupling 
schemes, involving also adjacent d spins, contribute to 
the MC wavefunction too. Such d — p — d interactions 
in the presence of the oxygen hole induce actually fer- 
romagnetic spin correlations between the d x 2_ y i electron 
at the "central" Cu c site and the nearest-neighbor d x 2_ y 2 
electrons, as illustrated by the Mulliken spin populations 
listed in Table I. 

The dressing of an O hole by a two-dimensional fer- 
romagnetic spin polarization cloud was previously sug- 
gested by Hizhnyakov and Sigmund [25] on the basis of 
Hubbard-like pd model Hamiltonian calculations. Never- 
theless, the dominant interactions are those between the 
d and p electron holes on a single CUO4 plaquette such 
that this object may be viewed as a ZR-like quasiparti- 
cle. As expected, the data from Table I also show that 
farther d—d interactions, between first and second order 
Cu neighbors to the central plaquette are antiferromag- 
nctic. 



III. RENORMALIZATION OF THE HOPPING 
INTEGRALS BY SPIN INTERACTIONS 

The effect of spin interactions on the nearest-neighbor 
and next-nearest-neighbor hopping integrals can be eval- 
uated by calculations on clusters of different sizes, where 
the O-mediated d—d magnetic couplings are included ei- 
ther for all or only part of the transition metal neighbors. 
For clusters that are large enough, we can take into ac- 
count both the ferromagnetic correlations between the d 
electron on the ZR plaquette and its nearest neighbors 
and the effect of the background antiferromagnetic cou- 
plings. 

The clusters we employ for our investigation always in- 
clude a 2-plaquette "active" region with one oxygen 2p 
hole plus a variable number of neighboring (undoped) 



plaquettes. The copper and oxygen ions on the two "ac- 
tive" plaquettes are modeled by all-electron split-valence 
basis sets of triple-zeta quality [Hj]. Depending on the 
size of the cluster, for the other cluster ions we either 
apply the same all-electron basis sets or use effective po- 
tentials for the core electrons. As mentioned above, point 
charges and total ion potentials are used for representing 
the rest of the crystal. For each cluster, we first determine 
the two equivalent CASSCF solutions where the 2p hole 
is localized either on the "left" or on the "right" plaque- 
tte of the active region. Since these states are obtained 
by separate SCF calculations, they are non-orthogonal 
and interacting. Non-interacting, orthogonal eigenstates 
can be obtained by State Interaction (SI, or CASSI) cal- 
culations [3(J. The effective hopping integral is half of 
the energy separation between these CASSI eigenstates 
[3l| . In terms of matrix elements between the separately 
optimized CASSCF wavefunctions, it can be written as 
t cS = (Hij - SijHu)/(l - Sfj), where H i3 and de- 
note Hamiltonian and overlap matrix elements between 
TV-electron states i and j, respectively. 

For all calculations we use a minimal active space, 
with one active orbital for each d and p hole. The total 
spin multiplicity is doublet because we always consider 
an even number of plaquettes (or Cu sites) and a single 
doped hole. During the CASSCF optimization for a par- 
ticular 2p-hole quasilocalizcd state not only the electronic 
charge is allowed to relax but also the spin configuration 
of the nearby Cu d x 2_ y 2 electrons is let to readjust, see 
for example the discussion in Refs. [H, HH, [H, [35|. Our 
approach for estimating the effective hoppings resem- 
bles in fact an idea proposed earlier by Eder and Becker 
|3(| , of setting up an effective tight-binding Hamiltonian 
H c s = ^efl ( c \ c j + h.c), where c] and c, consist of a 
sum of products of a "bare" fermion operator and a num- 
ber of spin flip operators. We note that for each hopping 
process within a "rigid" Neel lattice the energy of the sys- 
tem increases by 3 J/2 (35[. This quantity, about 0.2 eV 
(the value of the antiferromagnetic coupling constant is 
120—140 meV [37}), is comparable with the bare nearest- 
neighbor hopping of 0.4-0.5 eV, see Refs. 0, i, 1, 0, S 
and the results discussed below. By considering also the 
spin polarization effect described in the previous section, 
it can be argued that the magnetic energy associated 
with the disordered spins created at each hopping (3f| 
has the same magnitude with the bare hopping integral, 
i.e. the hopping of the 2p hole and the spin relaxation 
in the immediate vicinity imply similar time scales. Our 
procedure of mapping the ab initio data onto an effective 
tight-binding Hamiltonian including also spin relaxation 
effects is thus justified. 

pd ZR-like solutions where the oxygen hole is equally 
distributed over the four ligands of a given plaque- 
tte could not be obtained for 2-plaquette and linear 4- 
plaquette clusters with no distortions. For such clusters, 
the O hole has always the largest weight onto one of 
the ligands connecting two Cu ions. The smallest lin- 
ear cluster where the CASSCF calculations converge to 
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FIG. 2: (a) Sketch of the 10-plaquette [C1110O29] cluster used 
for evaluating the nearest-neighbor hopping t. The ions in the 
"active" region are shown in grey. The 6-plaquette [CueOi7] 
cluster is obtained by removing the four plaquettes on the 
same axis with the Cu ions of the active region, (b) The 8- 
plaquette [CU8O24] atomic configuration used for evaluating 
f. 



symmetric pd ZR-like solutions is a 6-plaquette cluster. 
Results obtained by CASSCF and CASSI calculations on 
the 6-plaquette linear cluster are shown on the first line 
of Table II. All-electron basis sets were applied in this 
case for each ion of the cluster [28|. 

The minimal-CAS SI nearest-neighbor hopping for the 
6-plaquette linear cluster is 0.310 eV. This value de- 
creases dramatically, however, when the cluster is en- 
larged such that it includes the other four CUO4 pla- 
quettes which are first order neighbors of the two active 
plaquettes. This is a 10-plaquette cluster, schematically 
drawn in Fig. 2 (a). In order to reduce the computational 
effort we used again, as in the 9-plaquette square clus- 
ter, the ECPs of Seijo et at. [H[ and Bergner et al. [l6| 
for the copper and oxygen ions outside the active re- 
gion, with the following basis sets for the Cu 3p,3d,As,Ap 
and O 2s, 2p outer shells: Cu (9s6p6d) /[3s2p2d] and O 
(4s5p)/[2s3p\. The CASSCF/SI estimate of the nearest- 
neighbor hopping integral for the 10-plaquette cluster is 
given on the last line of Table II. This value, 0.135 eV, is 
less than half of the estimate for the 6-plaquette linear 
cluster. As mentioned above, both charge and spin re- 
laxation are allowed during the SCF optimization. The 
Mulliken spin populations for the CASSCF solution with 
the O hole localized on the left-hand plaquette of the ac- 
tive region are displayed in Fig. 3. It can be easily noticed 
that the symmetry equivalent CASSCF ZR-like solution 
with the O hole on the right-hand plaquette implies some 
rearengements for the spin configuration at the neighbor- 
ing Cu sites. 

CASSCF and CASSI calculations were also carried 
out on a 6-plaquette cluster obtained by removing from 
the 10-plaquette cluster those four plaquettes on the 
same axis with the copper ions of the active region, 
see Fig. 2(a). The nearest-neighbor hopping matrix el- 



ement is approximately 0.17 eV in this case, see Ta- 
ble II. This value is again much lower than the num- 
ber extracted from the linear cluster. Still, since spin 
interactions are taken into account for only part of the 
nearest-neighbor plaquettes around the ZR-like quasipar- 
ticle, it is larger than the estimate for the 10-plaquette 
cluster. It is worthwhile noting that the value of this 
parameter changes only little when using effective core 
potentials instead of all-electron basis sets for the ions 
outside the active region, see the results for the [CU6O17] 
cluster in Table II. We also point out that the occupa- 
tion numbers of the cr p-d B and AB combinations on 
the ZR plaquettes and the Mulliken charge populations 
of the 2p atomic orbitals accomodating the oxygen hole 
are very similar for all the Cu-0 clusters discussed here. 
For comparison, Mulliken charge populations are listed 
in both Table I and Table II. Due to the symmetry prop- 
erties of the atomic configurations given in Table II and 
Fig. 2 (a), the p y ZR-like orbitals on a CUO4 plaquette are 
not equivalent, and we denote these orbitals as p v and p' y , 
where p y connects the two Cu sites of the active region. 

A different numerical experiment that confirms the ef- 
fects illustrated in Table II is to replace the eight Cu 2+ 
3d cations surrounding the two active plaquettes in 
the 10-plaquette fluster, see Fig. 2(a), with closed-shell 
Zn 2+ 3d 10 ions [H|. The CASSI estimate for t in this 
[Cu2Zn 8 C>29] cluster is 0.450 eV, larger than all values 
listed in Table II. For a [Cu@Zn4029] atomic configura- 
tion where only those four Cu neighbors are replaced by 
Zn 2+ 3d 10 ions which are not situated along the same 
line with the "active" Cu sites, t is lowered to 0.330 eV, 
close to the estimate in the linear [CugOig] cluster from 
Table II. This latter value is about two times larger than 
in the [CU10O29] cluster and indicates that the effect of 




FIG. 3: CASSCF Mulliken d-spin populations for a 10- 
plaquette cluster with one O hole. The ZR-like plaquette 
is represented with a filled square. The MSPs for the Cu 
and each O ion on this plaquette are 0.07 and —0.01, respec- 
tively. The equivalent CASSCF state used to calculate the 
effective nearest-neighbor hopping can be obtained by a re- 
flection through one of the symmetry planes of the cluster. 
For the two CASSCF ZR-like states, different spin popula- 
tions are associated with the nearby Cu sites. 
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interactions with and among Cu d spins which are not lo- 
cated along the direction of propagation is stronger than 
the effect of the collinear spin interactions. 

The next-nearest-neighbor hopping, usually denoted as 

was evaluated by calculations on an 8-plaquette cluster 
like that sketched in Fig. 2(b). We applied the same ECP 
basis sets as for the 10-plaquette cluster for the O and 
Cu ions surrounding the 2-plaquette active region, see 
above. All-electron basis sets were employed for the lig- 
ands in the active region [28|, but the two active cations 
were also modeled with ECPs [lij . The following ba- 
sis sets were applied for the 3p,3d,4s,4p shells at these 
transition metal sites: Cu (9s6p6d) /[3s3p3d] 15]. Our 
CASSI estimate for t' is 0.015 eV, which gives a t'/t ra- 
tio of 0.11 in the weakly doped La 2 Cu04 system. A ratio 
t'/t = 0.12 was obtained by fitting the ARPES data for 
an overdoped La2_ a; Sj , a ;Cu04 sample [39| and t'/t = 0.24 
for Bi 2 Sr 2 CaCu20 8 [l|. 

A few words are in place concerning the comparison 
between our present results and the results of other ab 
initio wavefunction-based investigations for copper oxide 
superconductors. A first issue is related to the nature of 
the active orbital space employed for the multiconfigu- 
ration calculations. For the 2-plaquette "active" cluster 
region we always consider a CAS with three electrons 
and three orbitals. These three orbitals are a bonding 
p—d linear combination, an essentially nonbonding d (or 
d—d) component, and an antibonding p — d orbital. The 
bonding p — d combination is usually [H, Q included in 
the inactive, doubly occupied orbital set, which gives an 
active space with one electron and two orbitals if other d 
orbitals around the 2-plaquette active cluster region are 
omitted. However, the occupation numbers associated 
with the pair of p — d bonding and antibonding orbitals 
in the 3-orbital CAS, about 1.8 and 0.2, respectively (see 
the discussion above and Fig.l(b-c)), indicate significant 
non-dynamical correlation effects [ill ] . As a zeroth-order 
approximation, the 3-orbital CAS seems to be then more 
appropriate for studying the character and the properties 
of a doped hole. Multi-reference, difference-dedicated 
configuration interaction (DDCI) [4Cj calculations with 
reference wavefunctions having either two or three or- 
bitals in the active space were performed on a hole doped 
2-plaquette cluster by Calzado and Malrieu They 
found that in the 2-plaquette cluster the DDCI estimate 
for the nearest-neighbor hopping is in fact quite insensi- 
tive to the size of the active orbital space. In addition, 
an iterative DDCI (IDDCI) scheme [41j was applied in 
Ref. for generating average molecular orbitals and to 
obtain an unbiased representation of the two states used 
to extract t. It was concluded that the description of the 
oxygen hole in the IDDCI wavefunction is only little af- 
fected when decreasing the size of the active space from 
three to two orbitals. 

For a 2-p laquette [CU2O7] cluster and all-electron ba- 
sis sets [28[ our CASSI estimate for the nearest-neighbor 
hopping integral is 0.650 eV, rather similar to the IDDCI 
value reported in Ref. 0, t = 0.575 eV [HI]. As already 



mentioned above, the "left" and "right" -localized O hole 
states are asymmetric in the 2-plaquette cluster. Using 
the same notations as in Table II, the Mulliken charge 
populations of the p y , p' y , and p x orbitals are 1.32, 1.74, 
and 1.68, respectively. 

A straight comparison between our CASSCF/SI re- 
sults and the DDCI/IDDCI data of Calzado and Mal- 
rieu is not possible. It is known 0, H^, [44j, however, 
that DDCI/IDDCI brings rather small corrections to the 
CASCI/CASSCF estimates for the nearest-neighbor hop- 
pings on 2-plaquette clusters in various cuprates, not 
larger than 15% [H, EH- This suggests that dynamical 
electron correlation is not essential in determining the 
magnitude of the hopping integrals, as already pointed 
out in Refs. [43l . |44| , for example. In any case, the spin 
polarization and relaxation effects at adjacent Cu sites 
evidenced by our calculations on clusters of different sizes 
are mainly related to non-dynamical electron correlation, 
and induce changes of the hopping parameter about an 
order of magnitude larger than the corrections brought 
by DDCI in 2-plaquette clusters. We also note that with 
the CASSCF/SI approach we are able to treat explic- 
itly charge relaxation and polarization due to the cre- 
ation of an oxygen hole on a given CUO4 plaquette. In 
the context of ab initio calculations for estimating ef- 
fective electronic-structure parameters such effects are 
sometimes referred to as dynamical repolarization effects. 
That a compact, transparent, and rather accurate de- 
scription is obtained when the most important configu- 
rations contributing to the many-electron wavefunction 
are expressed in terms of individually optimized orbital 
sets was previously shown for the calculation of magnetic 
exchange coupling constants [45| and core-level photoion- 
ized states [46j in transition metal oxides. 

An effective nearest-neighbor hopping parameter t can 
also be obtained by mapping the low-energy electronic 
states of the one-band Hubbard model onto those of a 
pd model. The pd Hamiltonian is given by the following 
expression: 

(ij)" 

+ ( d d L d ia + e P P\<yP3° 

+ U d 'Y^d\^d i] d\ i d ii + U p ^p] ] p n p] i p n . 

i 3 

t p d is here the hopping matrix element between p and 
d orbitals at nearest-neighbor ligand and metal sites and 
the first sum is over all pairs of such nearest neighbors. It 
is assumed that t p d takes the same value for all these Cu- 
O pairs. The second sum in the expression above is over 
all pairs of nearest-neighbor anions and t pp is the corre- 
sponding matrix element. We have chosen the phases of 
the orbitals such that the sign of each of the two hopping 
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matrix elements is constant. The difference between the 
on-site d and p orbital energies is in the electron repre- 
sentation A p d = e d — e p > and Ud > 0, U p > are the 
on-site Coulomb repulsion energies for Cu and O, respec- 
tively. For the parameters of our pd model we use typical 
values such as t p d — 1.5, t pp — —0.8, A p d = 3, Ud = 8, 
and U p = 4, in units of eV. The effective one-band Hub- 
bard model is given by 

+ U J2 c h c ^ c li c il ■ 

i 

t corresponds here to the effective hopping of the ZR-like 
quasiparticle between nearest-neighbor Cu04 plaquettes. 
U > is the "on-site" Coulomb interaction and we choose 
U=5 eV. 

Ground-state and excited-state calculations for the pd 
and the one-band Hubbard models are performed by us- 
ing DMRG techniques. For hole-doped clusters of 6 x 1 
and 6x3 CUO4 plaquettes we found that the lowest-lying 
states in the two models have the same quantum numbers 
and the same character. We can make thus between the 
two models a simple one-to-one correspondence for sev- 
eral low-energy states. The effective hopping parameter 
t is then determined such that it minimizes a function 

4 
n=l 

where AEff = Eff - Eg* (AE^ uh = E* uh - E$ uh ) are 
relative energies of the lowest four excited eigenstates in 
the pd (Hubbard) model. The ground-state is denoted 
with n = 0. For 6x1 and 6x3 clusters, each doped with 
two holes, our mapping procedure gives effective t values 
of 0.51 and 0.28 eV, respectively. This reduction of t is 
due not only to the inclusion of antiferromagnetic cor- 
relations on the adjacent chains but also to an increase 
of the nearest-neighbor spin correlations along the chain. 
Within the pd model, the nearest-neighbor spin correla- 
tion function changes from —0.103 in the 6x1 chain to 
—0.152 for the central chain in the 6x3 ladder. For more 
chains this quantity is only slightly modified, to —0.158 
for example in a 6x5 cluster. The model-Hamiltonian cal- 
culations confirm thus the renormalization effect found 
by ab initio quantum chemical methods. If six holes are 
doped into the 6x3 cluster, i.e. the hole concentration is 
kept constant with respect to the 6x1 2-hole configura- 
tion, the effective hopping is 0.29 eV, which shows that 
for this range of hole concentrations the doping depen- 
dence is very small. 



IV. CONCLUSIONS 



The interplay among charge, spin, and lattice degrees 
of freedom and the formation of composite spin-lattice 
polarons were studied in cuprates (and other transi- 
tion metal oxide compounds) by many authors. In lay- 
ered copper oxides, the spin mechanism responsible for 
polaronic-like behavior is associated with the fact that a 
charge carrier hopping to an adjacent site must disrupt, 
at least at relatively low doping, the antiferromagnetic 
spin configuration in its vicinity. Such effects were inves- 
tigated by calculations based on one-band Hubbard and 
t-J models, see for example Refs. [13, SI, IH, US HH, H3j • 
These studies show indeed that the background antifer- 
romagnetic correlations lead to band narrowing and a 
strong enhancement of the electron-phonon couplings. 
We are able to provide here ab initio numerical estimates 
for the effect of such interactions on the effective hop- 
ping integrals. The results can explain the relatively low 
values deduced for the quasiparticle hopping parameters 
from fits of the photoemission data. 

Our ab initio calculations also show that an oxy- 
gen hole induces ferromagnetic correlations between the 
d electron on the ZR-like plaquette and the nearest- 
neighbor Cu d spins. Such effects are missing in one-band 
Hubbard or t— J models. At this moment we are not able 
to quantify the strength of these ferromagnetic interac- 
tions in the form of effective magnetic coupling constants 
or spin correlation functions, at least not at the ab initio 
level. Nevertheless, this spin polarization effect is one of 
the factors that determines the characteristics of the hole 
motion. It might also offer an explanation for the disap- 
pearance of long-range antiferromagnetic order at doping 
levels which are much lower as compared to the elec- 
tron doped cuprates. It is interesting that for an isolated 
chain of plaquettes, "local" ferromagnetic correlations in 
the direction of propagation would tend to facilitate the 
motion of the hole, although the fact that the carrier per- 
turbs the magnetic couplings beyond nearest neighbors 
results in a stronger and opposite effect. In addition to 
that, the effective hopping parameter depends strongly 
on interactions involving spins on adjacent chains, i.e. 
Cu d spins which are not situated along the direction of 
propagation of the hole. Qualitatively, this is also con- 
firmed by mapping the solutions of a pd model Hamilto- 
nian onto those of an effective one-band Hubbard model 
for clusters consisting of chains or ladders of plaquettes. 
The inclusion of the adjacent antiferromagnetic chains 
in a three-leg ladder determines a reduction of the effec- 
tive nearest-neighbor hopping of about fifty percent in 
comparison to the case of a single chain. 
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